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1. Fourier lias shown that if any function f(t) (=.y say) of a 
variable t be such that 

/(/)=/(/+ r )=/(/+2r) = etc., 

where r is a constant, that is, if /(/) be periodic in t , of period r, 
then /(/) can be expressed as the sum of a constant and a series 
of terms called harmonics, each of the form 
a p sinp(u>t - 6 p ), 

where p has the values 1, 2, 3, 4, etc., 
and to = 27r/r. 

The number / is called the order of the harmonic, a p its ampli¬ 
tude, and Op its phase. 

If, in addition,/(/) be such that 

_ /(')= -/('+r/2), 

then it is easy to see, by substituting /-fr/2 for /, i.e., to/+ 7r for 
to/ in 

y=a o + ’2a p &inp(<at—0p) i 
that in order for y t to be = —yt+rj 2 
a {) = 0, = a± — 0, etc. 

Hence in this case the constant term vanishes and the har¬ 
monics, of which /(/) is the sum, are all of odd order. When 
such is the case /(/) is called an odd periodic function. This is 
the type generally met with in alternating electric current in¬ 
vestigations. 


1 Appendix to the paper: “Preliminary Account of a Wave Tracer and Analyzer.” 
Phil. Mag., Nov., 1903. 
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2. If we define the nth component (C n say) of a periodic 
function /(/) of period r as the periodic function which is the 
sum of those harmonics of f(t) whose orders are n , 3;/, 5 n, 7//, 
etc., tlien 

2„o,=/»-/(/+£)+/(/+*£)- .... 

-/(/+-2—1^). (I.) 

For if we represent the expression on the right of the above 
equation by \we find by substituting successively for /, 
/ + r/2« and t-t-rjn in it, that 

«'>=- K '+§,)=*('+;)■ 

Hence i^(/) is an odd periodic function of period r/;/, that is to 
say, if 

/(0 = <*o + 2tf p sin/(w/- 0 p ), 
where p— 1, 2, 3, 4, etc., 
then i^(/) is of the form 

i p(t) = 'Zbqsmqn(u)t - f3 q ), 
where q— 1, 3, 5, 7, 9, etc. 

In evaluating ij/(t) therefore, only those harmonics whose 
arguments are mat, 3/W, 5 mat, etc., need be considered. Ne¬ 
glecting all other harmonics in the different /functions that make 
up ij/(t), we find that the remainders in the 2 n terms 

./('). -/(' + 2 ^)- ^(' + 2 ^)> etc ‘> 
are all equal, and that each remainder is the «th component of 
/(/), lienee 

'l'(t) = 2«C n . 

3. If /(/) itself contain only odd harmonics as in the case of 
alternate current periodic functions, then 

/«=-/('+;> 

and equation I., §2, reduces to 

=/(/)-/^+ 2^)+ • • • +/(*' + 71 ~ *2^)- 

The operation on /(/) mathematically represented on the right 
hand side of equations I. or II., is practically performed on 
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alternate current waves by the wave tracer and analyzer 1 
designed by the author. In the simplest case, when n = 1, the 
wave tracer gives the first component of the periodic quantity 
operated on, which in the case of alternating electric currents is 
the full wave. By the movement of two pairs of brushes n can 
be made 3, or 5, or 7, in which cases the analyzer will give the 
3rd, 5th, or 7th components of the wave respectively. 

Now, in practical investigations with this apparatus on alter¬ 
nating current waves whose harmonic expressions were required, 
it was found much better to obtain by its means only the full 
wave trace, and then by an arithmetical process identical with 
the action of the analyzer and indicated by equation II. above, 
to obtain the 3rd and higher components of the wave, and thence 
to deduce its harmonics. 

This method of harmonic analysis was drawn attention to in 
the paper already quoted, and though based on a different 
formula to that of Wed more, 2 is practically similar to his. It is 
more suitable, however, for waves containing only odd harmonics, 
and as I have had considerable experience in its use during the 
last two years and have found it both expeditious and accurate, 
it is possible that a short account may be of value to those 
interested in alternating current work. 

4, In wave graphs it is more convenient to use angular abs¬ 
cissae x where 

x = uit=2irtjr. 

Making this substitution in the equation y = /(/) it becomes 
y = g(x) say, where 

A’d) =£•(*+2 tt), 

and if f(t) is an odd periodic function as in the case of alternate 
current waves which we are now considering, 

= - g(x + tt) =g(x H- 2 *). 

Substituting g(x) for/(/) in equation II, it becomes 
>‘C n = g(x) -°(x + +g(x + 2 njn) - . . . 

+g(x + n - 1 tt\h), 
from which we conclude that, if 

1 Lyle: “Preliminary Account of a Wave Tracer and Analyzer.” Phil. Mag., Nov. 
1903. 

2 W'edmore: Journal Inst. Elect. Engineers, vol. xxv., p. 224 (1S96). 
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y 0i Ji> J* 2 >- Jn-i be n equi-spaced ordinates that exactly in¬ 

clude half the wave, i.e., ordinates corresponding to the abscissae 
.r, x + ir/n, ,v + 2 7 r/«, . . . . x + h-Itt/h respectively, and 

called e.s. ordinates in the sequel; and if N 0 , N 2 -N n _i, 

be the ordinates of the wtlv component C n whose abscissae are 
the same as those of j' 0 , y v y 2 ,- - -y n -1 respectively, then 

• • +^ r n-i = «N 0 = =«N„_i 

when n is an odd number, and 

J0-J1+V2- • • • -J'n -1 = 0 

when 11 is an even number, as we are now considering odd 
periodic functions only. 

Thus from n e.s. ordinates of the original half wave we obtain 
only one ordinate per half wave of C n , so that in order to obtain 
111 e.s. ordinates per half wave of C n it is necessary to have inn 
e.s. ordinates of the original half wave. 

For instance, to obtain 3 e.s. ordinates of C n we must measure 
3 11 e.s. ordinates of g(x). Let these be 

J'oj Jl) J'l J'b>‘ ” J'3n-1} 
and let the corresponding ordinates of C n be 

L 0) 1) -^2) -^3? 

then 

Jo-Vi+ye- • • +^8u-3=«N 0 = —«N 8 =;iN b = = //N 3n _ 3 

)\ — y\-\-y \— • • +J’3h- 2 — — ;/N 4 — 7/N 7 =: = «N 3n 2 

J’z-y^+y^- • • +J , 3n-i = »N a =—»N 5 =»N fl = =;/N 3n _i 

Subtracting now the ordinates of O u so obtained from the cor¬ 
responding y ordinates, we obtain a new set of ?>n e.s. ordinates 
which are those of the original half wave with its ;/th component 
removed. 

5. In practice it will generally be sufficient to determine the 
1st, 3rd, 5th, 7th and 9th harmonics (H x H 3 H- FL say). 
This can be done with considerable accuracy when 15 e.s. ordin¬ 
ates of the original half wave are given. 

Thus if these be 

y» )\y y& ♦ • • j’u 
corresponding to the angular abscissae 

*^'0J • * * ^14 

where x x - x Q = x 2 - x 1 == . . . = x u - x ls =tt/1 5, 

and if z 0) g l9 z.,, z S} z 4) be 5 e.s. ordinates of the half wave of 
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C 3 , then 

3*o =J>' 0 ~y 6 + Jio = - 3z 5 = 3z 10 
3*i —)’i ^ ~ 35 6 =3s u 

3*4 ~ + J'i4 = ~ 3 s 9 = 3s 14 , 

and if ;/ 0 , 7/ l5 7/. 2 be 3 e.s. ordinates of the half wave of C 5 , then 
0Uq — j ’ 0 y z y$ y g -f- y^— ^^3 = 5 77,5 = 5 u 9 — 5 u u 

5 "x=J’i— J 4 +/ 7 — J’io+Ji3= -5u 4 =5u 7 = -5u 10 = 5u u 
= y* —y s + >8 - J 11 +7i4 =-5u 5 = 5u 8 =- 5u n = 5* I4 
the figure subscribed to each ordinate indicating the abscissa to 
which it corresponds. 

Now the full wave 

C x = H x + H s -f- H 5 -j- H 7 -f- H 9 -j- etc. 
and C 3 = H 3 d-H 9 4 -H 15 

C 5 = H 5 + H 16 , 

so that if H 13 be neglected, and the sums of the corresponding 
ordinates of C 3 and C 5 be subtracted fram those of C lf the fifteen 
remainders are ordinates of 

H 1 + H 7 + 

i.c.y of H x , if we neglect H 7 . 

If H 15 cannot be neglected it can at once be removed from C 5 
before subtracting from C 1? for as it is (q.p.) the 3rd component of 
C 6 , of which we have 3 e.s. ordinates 77 0 , 7 / 2 , its three corres¬ 

ponding ordinates are ? 0 , —/ 0 , / 0 where 3i 0 =?t 0 — t/ 1 + 
hence H 5 will de completely given by 
c 0 , <r 1? c . 2 where 

7/ 0 —/ 0 , Ci = 77i + /(j, ^*2 — 7/ 2 — 7 0 . 

H 16 can now be taken from C 3 , thus 

*0 “ *1 d" ^05 *2 __ ^0) *3 d" *4 — *0 

are the 5 e.s. ordinates of H 3 +H 9 . 

In order to determine H 3 and H 9 it will now be necessary to 
plot the 5 ordinates of H 3 -f H 9 , measure off 6 e.s ordinates from 
the smooth curve drawn through them, and from these determine 
their first component, that is 2 e.s. ordinates of H 9 . These will 
completely determine H 9 if H 27 etc., be neglected, and by sub¬ 
tracting them from the corresponding ordinates of H 3 -+-H 9 6 e.s. 
ordinates of H 3 are obtained. 

If H 7 cannot be neglected it will be necessary (if the orig¬ 
inal wave trace is not available) to plot the 15 ordinates of 
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Hj + H 7 obtained above, and from the smooth curve drawn 
through them to measure off 14 e.s. ordinates. From these, 2 
e.s. ordinates of the half wave of H 7 , and which determine H 7 , 
can be obtained. By subtracting these from the corresponding 
ones of 14 corrected ordinates of H x are obtained. 

G. It now remains to determine the amplitudes and phases of 
the harmonics of Cj from their ordinates which we have obtained. 
It is easy to show that 

-|sin 2 0 -f sin 2 ^ + sin 2 ^0 +—J+ • * 

+ sin 2 (0 + u - 1 J)} = 1, 

from which we conclude that the square root of twice the mean 
of the squares of n e.s. ordinates of half a sine wave is equal to 
its amplitude. 

Hence, with the help of a table of squares or of the quarter 
squares given in most sets of tables the amplitudes of H l5 H 3 , 
etc., can be quickly determined. 

[The rule that the amplitude is equal to 7 r /2 x mean of the 
ordinates is only sufficiently accurate when a large number of 
ordinates is taken.] 

If « 0 , a 2} . . . a u be the ordinates we have found for 

Hj corresponding to the angular abscissae x 0 , x l9 x 2 . . . x u 

respectively, and if /q, a be the amplitude and phase of or in 
other words, if 

H! —/qsi n (cut — a), 
then any of the equations 
sin (x 0 -a) = ajh 1 
sin(jCj - a) = a 1 /h 1 
sin (^2 - a) = a 2 //i 1 etc., 

would determine a, provided the ordinates « 0 , <q, a. 2i etc., are 
exactly those of a sine wave. 

In practice, however, small upper harmonics will invariably be 
left in a 0 , a 2 , etc. [it may not have been thought worth while 
to remove H 7 ], and though their amplitudes may be negligably 
small, yet they might cause considerable error in the value of a 
when determined from only one of the above equations. Hence it 
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is advisable to obtain four values of a from the first four ordin¬ 
ates on the rising side of the wave and four from the last four 
ordinates on the falling side, and take the mean of the eight. In 
this way we can to a great extent eliminate any error that might 
arise due to a harmonic even as low as the seventh not having 
been removed. 

In a similar way the phases of H 3 , H r , etc., can be determined, 
but it must be remembered that if, for instance, 

H 3 = /z a sin3(oj/ - /?), 

and if /; 0 , b x , . . b 4 are the ordinates of H s corresponding to 

the abscissae x Q , x l4 . . x 4i then 

sin3(.r Q - fi) = b 0 //i 3 etc. 
similarly, if 

H 5 =// 6 sin5(a>/ - y) 
with ordinates 4 , 4 , 4 , 

then sin5(.T 0 - y ) = c 0 /h 5 . 

7. The wave to be analyzed may be given in either of two 
ways. We may have the complete trace of it obtained by the 
author’s wave tracer by the photographic method, or by any 
form of oscillograph that gives a trace of the wave form ; or we 
may have the values of a definite number only of ordinates per 
half wave, such as would be obtained by the author’s wave 
tracer by the galvanometer and scale method. 

From the wave trace the complete harmonic expression can 
theoretically be obtained, but the impossibility of accurately 
measuring on the photograph, without elaborate apparatus, the 
different ordinates required leads to great inaccuracy in the 
result. 

From a given number of e.s. ordinates only an approximate 
analysis can be obtained, more approximate, of course, as the 
number of ordinates is greater. When, however, each individual 
ordinate has been obtained with the accuracy of which the galva¬ 
nometer and scale method is susceptible, the analysis obtained 
from fifteen such ordinates is much more reliable, as far as the 
harmonics up to the 9th are concerned, than that determined 
from any photographic trace. 

I will therefore illustrate the method by applying it in full 
detail to the analysis of the wave whose 15 e.s. ordinates are 
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given in row 5 of Table I. Every figure necessary in the calcu¬ 
lation will be given. 

The first row of figures in Table I. are the abscissae .v 0 , x lt etc. 
to which the given ordinates correspond. Space for three rows 
of figures is left, and then the 15 given ordinates are written 
down. These are divided into three sets of five each, and the 
numbers of the middle set are subtracted in order from the sums 
of first and last set, giving five numbers which are the corres_ 
ponding ordidates of 3C S . Space for two or more rows is left 
and the given ordinates are now written down as in the table, in 
two rows of six each and one row of three, in order. The 
columns formed are added and the last three of the sums are 
subtracted from the first three, giving three ordinates of 5C 5 . 
The first of these minus the second, plus the third, gives one 
ordinate of 1 5C* whose other ordinates are got by alternating 
the sign. Subtracting 5C 15 from 5C 5 we obtain 5H-. Having 
obtained C 15 we now subtract 3C 15 from 3C 3 and obtain 3(H 3 -b H y ). 

Above the given ordinates write those of C H with signs changed 
(row 4), and above these write those of H. with signs changed 
(row 3). Add rows 3, 4 and 5 to get row 2, in which are the 
ordinates of Hj + H 7 H^ etc. Neglecting Il 7 , H 1]? etc., as is done 
the analysis in Table I., we may consider the figures in row 2 as 
the ordinates of H,, and neglecting H y we may consider the 
figures in row 11 as the ordinates of 3H 3 . 

The first 15 numbers under Amp. H 1 are the quarter squares 
of the ordinates of H,. Twice the sum of these is divided by 15, 
the number of ordinates, and the quotient is found to be the 
quarter square of 987. Hence /q, the amplitude of Hj, is 987. 
Similarly for the amplitudes of H 3 and H r> . 

Under the heading “phase of H 1} ” in the first column under 
sines, are the quotients got by dividing the first four ordinates 
on the rising side of Hj and the last four on the falling side of H : 
by /q; in the second column under angles are the corresponding 
angles, and in the third column are the eight values of 12 ° — a 
deduced. The mean of these 2° 2' when subtracted from 12° 
gives the crossing point or phase of H x as 9° 58'. Similarly for 
the phases of H 3 and H 5 . It will be noticed that at the crossing 
point determined for H s , H a crosses down, which is expressed 
analytically by writing its amplitude negative. 
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8. It will be noticed in the determination of the phase of Hj 
in Table I., that the eight values of 12° — a differ considerably 
from each other, indicating the presence in what we there take 
for Hj of a considerable upper harmonic, probably H 7 . Jn order 
to determine H 7 fourteen e.s. ordinates of the half wave are 
required. Tf the wave trace were given these could be measured 
off from it, but if, as in the case we are considering, only 15 
original ordinates are given, it is necessary to plot the 15 ordin¬ 
ates of Hj-j-H 7 obtained in Table I., and from the smooth curve 
drawn through them to measure off 14 e.s. ordinates. This has 
been done and the values obtained are given in row 4, Table 
II., as well as the calculation necessary for the determination of 
H 7 and its elimination from 

What is called the amplitude of H x in Table I. is really 
v /2 R.M.S. (H 1 + H 7 ), To get amp. H, it is better to remove 
the effect of H 7 by treating it as a correction, thus avoiding 
error that might be introduced in the plotting. This is easily 
done, since 

M.S(H 1 + H 7 )=^-V/ 

hence h x (corrected) = ' Amp.(IIj-f-H 7 ) 2 —// 7 2 
— V h* (uncorrected) — // 7 ‘ 2 . 

In Table II. the corrected crossing point of Hj is determined, 
and it is seen to differ in phase only by 2 minutes from the value 
obtained in Table I. 

The difference between the four values of 3(24°— j3) when 
determining the crossing point of H 3 in Table I. point to the pre¬ 
sence of a ninth harmonic, which exists as a third component in 
II 3 4-H 9 . H 9 can, if desired, be determined by plotting the 5 
ordinates obtained in Table I., measuring off from the curve six 
e.s. ordinates and proceeding as before. It will be found that 
H 9 =3sin9(a>/—13°). 

9. In Table III. is given most of the work required for the 
determination of the first six harmonics of a complete wave that 
contains harmonics both of odd and even orders. Twenty-four 
e.s. ordinates of the full wave are taken. This number is speci¬ 
ally suitable, as it enables us to determine directly Cj, C 2 , C 3 , C 4 
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and C 6 . To determine C 6 , replotting will have to be resorted to 
if the full wave trace be not available. 

At the top of Table III. are written the 24 given ordinates 
under their corresponding abscissae. From these ordinates the 
constant term of f ( t) has been removed. This can be done by 
aid of the formula 

/(0+/(Hr/«)+/(/+2r/«)+ . . . +J{t + n- 1 T /n) 

= n\a 0 + a n sin//(a>/ — 0 n ) + ao n s\n'2n(wt—0z n ) 

-f rt 3M sin3«(aj/— 0 3n ) -f etc.] (TII.) 

which can be easily established by the method used in § 2. 

From this formula we see that the mean of n e.s. ordinates 
embracing one period of a periodic function is equal to its con¬ 
stant term, if its «th, 2//th, etc., harmonics are neglected. 

Returning to Table III., we add the second twelve ordinates 
with their signs changed to the first twelve, in order, and obtain 
12 e.s. ordinates of 20^ i.e., of 2[H 1 + H 8 + H 5 + ]. (See equa¬ 

tion I., § 2). 

Subtracting these from twice the given ordinates, those of 
2[H 2 + H 4 -1- H 6 + ] are left, and the remainder of the work 

proceeds as in Table I. 

2[H 2 + H 4 + H r> +- etc.] could be obtained directly from the 24 
given ordinates by adding the second 12 to the first 12 of them, 
in order. (See formula III., § 9). 

The amplitudes and phases of the different harmonics were 
determined as in Table I., but the figures necessary in their 
calculation are not given. 


The following are interesting applications of the above method 
to more general harmonic analysis. 

10. To obtain the harmonic expression for the odd periodic 
function whose graph for half a period is the sides of an isosceles 
triangle of altitude h. See Fig. 1. 

Taking o and i r as the abscissae of the extremities of the base, 
relative values of any number of e.s. ordinates can be written 
down, and any component at once obtained. Thus, 30 e.s. ordin¬ 
ates would be 0, 1, 2, 3, . . 14, 15, 14, . . 2, 1, and these 
correspond to an altitude 15. 
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It will be found that all the components (i.e. 3rd, 5th, etc., in 
this case) are the sides of isosceles triangles passing through the 
origin, and that the altitudes are 

3 2 , /z/5 2 , — /z/7 2 , etc. respectively. (See Fig. 1.). 

(The same can be quickly arrived at geometrically). 

Hence, if the full wave or Cj be represented by 

CJ X = <T x sin(to/ — Of) -t-tf 3 sin3(co/ — 0 3 ) -|-rt 5 sin5((D/— 0 5 ) + etc., 
its third component C 3 is 


= — Jo[^ 1 sin(3a)/—^ 1 )+z? 3 sin3(3a)/— 6 Z ) 4-rt 5 sin5(3a>/ — 0 5 ) 

o~ 


+ etc.], 
and its fifth component C 5 


1 

o 


) [z7 1 sin(5o)/ — 6i + Z7 3 sin3(5a)/— 0 3 ) + zz 5 sin5(5o)/— 0 5 ) + etc. 


] 


and so on, but by definition C a and C 5 are also given by 
C 3 =Z7 3 sin3(a)/— Of) + rt 9 sin9(u>/— Of) + 

C 5 = z7 5 sin5(co/— 0 r ) + « 15 sin 1 5(oj/— 0 15 ) + 
hence, identifying the expressions for the same components, we 
find that 


a x — — 3 2 tf 3 =5 2 rt 5 = — 7 2 zz 7 = etc., 
^=303=505=707=etc., 



